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A microscopic derivation of the Born—Huang relations 


between the atomic force constants 


By Lars T. HEpIN 


SUMMARY 


M. Born and K. Huang have derived a general relation [1, 2] between the elastic constants 
and the atomic force constants by considering long acoustic waves in an infinitely large 
lattice. The symmetry properties of the elastic constants lead to certain linear relations be- 
tween the atomic force constants that could not be predicted from a microscopic treatment. 

It is shown here that by taking the finite extension of the lattice into account the Born— 
Huang relations readily follow in the microscopic treatment. 


1. Introduction 


According to Hooke’s law the relation between the stress tensor, o,;, and the 
strain tensor, ¢,3, is linear, 


Op = > Cute (1.1) 


The elastic constants C,s,,, have the symmetry properties 
Cap, to CBee, uv Cuy, aps (1.2) 


The subgroup of the permutation group which leaves c,,,,, invariant will be 
called G and instead of (1.2) we write 


tiie @ Gs (1.3) 


There are in general 21 independent elastic constants. 

When the solid is under initial strain a small change in the stress tensor 
may still be linear in the changes of the strain tensor but the constants of 
proportionality, the “linearized elastic constants”, no longer belong to G (see 
appendix). 

The microscopic theory of elasticity is built on the concept of interatomic 
forces. These are derived from an atomic potential energy usually taken in the 
adiabatic approximation. The second derivatives of the potential are called force 
constants. Born and Huang have derived a relation between the force constants 
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and the elastic constants by comparing the microscopic equation for long waves 

with the equation for the acoustic waves as given by the theory of elasticity.t 

We conclude this introduction with a brief survey of the Born-Huang results 
and references are given to equations in their book [4]. 

‘ . l 

The atomic potential © is a function of the atomic coordinates 2, (,) where 

a is a cartesian index, J indicates a lattice cell and k a specific atom of the 

; l s : 

cell. They consider small changes of the coordinates from (;) a L, G, + & (k) 

=1 


r 
l 
corresponding to equilibrium positions, to @ (;) +% ( 4) The corresponding change 


in the potential is 
Wess oe (re) (;)« vt (1.4) 
“WEY Bee ONE ao ee ee ; 


where the quantities © are called force constants. As they consider an 


I’ 
PN bh! 
infinite lattice the absence of the first order terms in A@® does not exclude the 
possibility of a finite strain. On account of the lattice symmetry 


11! panty apy 
a ee | k k’ ) 


It is convenient to use the notations 


l 10 
ap (: v) = Pup 6 2) 
Ld ha 10 1 0 
d =F =f —fFr 


The method of “long waves” gives (compare eq. 28.7), 


Say, pat Sea, py = 2 [a B, yAl+ (ay, BA) + (ad, By). (1.5) 


Here Sz,,p, ls a coefficient in the elastic energy density (4.2) and we show in 
the appendix that it also is a linearized elastic constant. The symmetry pro- 
perties of the coefficients in (4.2) are 


Sap = Spas (1.6) 
Tope, (1.7) 
where sos (vis ie, pa— Oxp Savi (1.8) 


1 These relations may also be derived by the method of deformation fields [3]. Notice that 


th ticl ; : : i ; ; F 
3 ; fe a “ie by Begbie and Born contains some mistakes which were corrected in the article 
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The terms of the right member of equation (1.5) are (v,=volume of lattice 


cell), 
1 l l l 
laB, yA] = Bo, eye Po bed av, re) vy ee) (1.9) 


and («y, A) which is a complicated non-linear expression in the force constants. 
When the magnitude of the vector 1 increases the force constants D2 c Al 


rapidly tend to zero. The lattice sum (1.9) is expected to converge for all solids 
except the piezoelectric ones. 
The microscopic treatment of the infinite lattice shows that 


[xB, yA]=[Ba, yA]=[aB, Ay], (1.10) 
(«B, yA) EG. (1.11) 


Using (1.8) we may rewrite equation (1.5), 
[T'ay,aa— (cys B.A)\ + [Tar py — (HA, By)] = 2 {le B, y A] — bap Sya}- (1.12) 


The left member is invariant when we permute «yf into ya jf. Demanding 
this invariance of the right member we obtain new relations which together 
with eq. (1.10) may be writen 


Caf, ya € G, (1.13) 
where Cap,ya =[«B, yA] — bap Sya- (1.14) 


The tensor 
[a B, yA] = [dap Sya — bya Sap] + } [Bap Sya + bya Sap] + Capya 


has 26 independent components. This is easily seen since the first parenthesis 
on the right hand side may be described by 5 parameters and the remaining 
terms by 21 parameters. 

The relations (1.10) show that [«6,yA] has at most 36 independent com- 
ponents and thus the properties of the coefficients S,,,,, infer 36—26=10 new 
linear relations between the force constants. These are the “‘Born—Huang rela- 
tions’ (28.32)—(28.35). 


2. Invariance relations in a finite lattice 


We first derive some relations that are valid for an arbitrary system of N 
particles interacting via a potential (z,, @, ... Zy). We use a simpler notation 
than in the preceding section. A Greek index will specify the cartesian com- 
ponent and a Latin index the particle, for example; (#;)z=ig. The derivatives 
are written 


oO 
pha 2es Oy. =7—— > ote. (2.1) 
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We immediately have the relation 
Dijop = Djipe 
The potential is invariant under translations and rigid rotations, 


 (&,, Gp, ..., dy) =O (& +4, 44, ..., By +a) =O (PH, PH, ..., Py), 


where P is an orthogonal 3x3 matrix. 
Taking a partial derivate of (2.3) we obtain, 


Dig (Ey) «+19 Ey) = Dig (4, +4, ..., Eyt+ a) = 2, Py, (Py, See) eee 
We make a Taylor expansion of (2.3) and (2.4) in the components of @, 
O (Z,, .-2) = (4, ...) + 2, Pie de + ee 
Diz (Z,, ---) = Dig (By, -23) 4 2, Pires ap + a. 
Hence > oO, = 2: Oj;.3 = 0. 
When P corresponds to a small rotation it may be written 
Pp =Oap + Wapt...3 Wap= — Wpa- 
We make a Taylor expansion of (2.3) and (2. 4) in wags, 
® (€,, ...) = (&,, cas, Di_. Lig Wap t ++, 
Di. (#4, ...) = 2. (Opa + Opa + *+) (Din (%1,-.-) + 2 Dijuy Xjp Wyp + -*+) 
= Dj, (2, 06.) + ~ Dijay Lig Wyp + 2, Dip Oya + ++ 
Hence 2 O;,, mip = > Dig Lig, 
2 Dijay Lig + Diy Oxp = 2, Dijap Ljy + Dig Say 


Introducing the quantities 


Tap = 2 Di. Xi, 
v 


ny 
V 


1 
Geuvp ~ TF 2 Dijx8 Lip Lip, 
V Gj 


(2.2) 


(2.3) 


(2.4) 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


(2.10) 


where V is the volume of the convex polyhedron spanned by the particles, we 


obtain from (2.7) 
Tp = T Ba 
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and from (2.8) after multiplying by 2;, and summing over i 

Aeupy + Op T uy = Depp + Oey L pyu- (2.12) 

From (2.2) and the definition (2.10) of dy. we find 
Geyop = Fprpice- (2.13) 

The equations (2.12) and (2.13) may be combined to 
Das € G, (2.14) 
where Dayyp = Ceuyp — Ox6 T uv. (2.15) 


We now specialize our results to the case when the particles form a space 
lattice. In the previous notation the positions are 


From (2.2), (2.5), (2.9) and (2.10) we have 


0.2 € =) = Opa (;: nt (2.17) 
z(2)-z0e((4)-0 os 
ron} 30 (i)a() ew 
dauvp =F ap OD. c | Lp ( 4 My ( ), (2.20) 


Using (2.18) we obtain the identity 

ti Ll’ Ll’ el (;) (;) (;) (_.)} 
ita: ©, , vl, + ty rTf- 

pp P ty) ee bei) ie (ci) We ont wea Wma Vg seca 4) ca 
(2.21) 
The right member is equal to — V (dz + dup). The left member can be written 

Vel (’4+l’jeL i SEG Ms 1K 1H [’ . 
‘ 2.22 
> 2 2 a ( k ) al u)® (, x) oe 


The notations l’€Z and (l’+l’’) €L indicate that the vector components satisfy 
the conditions in (2.16). 
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kk’ kk’ 
(2.22) over 1’, & and k&’ gives —2u,[«B, pv») according to equation (1.9). In a 
finite lattice this only represents an approximation which becomes better the 
larger the distance between the point Z(l’) and the boundary of the lattice. 
Under certain! conditions we obtain 


l’ a8 te ih fee : ; 
If the lattice is infinite then .s( )= 0.0 and the summation in 


—2v,[«B, wry] N+ O(N) = —V [dayep + devyp)), (2.23) 


where N=N,N,N, is the number of lattice cells. Since we consider N as very 


large we write 
[oe B, wr) =$% (dayrp + Carus) (2.24) 


neglecting the term @(N~*), From (2.15) and (2.24) we obtain 
fapur = 3 (Daye ze DD egatls (2.25) 
where fu = (eB, WY] — bap Ty (2.26) 


Since Dawg EG and fupur = fea = fap, We Obtain from (2.25), (compare the deriva- 
tion of (1.13)), 
fap © G. (2.27) 


This shows that [«/, uv] obeys the Born—Huang relations. 
The essential steps in the proof are 


(i) rotational invariance of D gives Dy € G; 
(ii) neglecting certain terms associated with the surface gives a relation be- 
tween [«f8, wv] and Dwg, equation (2.25). 


3. Remarks 


The results in the previous section were obtained by considering a finite lat- 
tice whose dimensions became infinitely large instead of introducing the infinite 
lattice from the outset (compare the recent work by E. Stern [5]). It is con- 
venient to denote quantities referring to the limit case of an infinite lattice 
with a superscript °. 

We observe that in the derivation of the Born—Huang relations we only used 


| =0 but not 


: l : 
necessarily 2 ( ) =0, which represents a slight extension of the validity of the 


Born-Huang relations. 


the translational symmetry of the lattice. Then we have >, 02 


1 The difference between the force constants in a finite and an infinite lattice must de- 
crease sufficiently rapid with the distance from the boundary, and the lattice sum in (1.9) 
must converge sufficiently fast. The sufficiency conditions are however very weak. Thus if, 
for example, in (1.9) the sum over k, k’ and 1 satisfying |#(l)|<r gives a remainder term 
smaller than r~*, €>0 they are fulfilled. 
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We further note that the present technique makes it possible to obtain the 
relation between the elastic constants and the force constants by considering a 
homogeneous deformation instead of long waves. This has the twin virtue of al- 
lowing the complete derivation to be made with a finite lattice and being simpler 
than the method of long waves. We give a sketch of the derivation. 

Consider a small homogeneous deformation from the equilibrium positions 


(7)- 2, l,d,+2(k) to #(,)+a(;). 


where Ux (;) = > Uap Xp (;) +Ug(k), Usp > 0. (3.1) 
B 


The vectors u(k) are determined from the condition that 


Lie if 
0 = 
ee f v) oe (:) Y 


The change in the potential ® becomes, 


l l f 
=e (;) 2»(;) m+ 3% (;)me 


Pave seth i: 11 () : 
i k 
mM ee Pep (; i) os (;) ‘ i) Sea WA Pug (; ki’) "#\ 4p) Me up (hk) 
aBuy 


au 
tidy Don Dap ( ) Ua (k) Up ( .) . (3.2) 
wR kk k 


The second and the last terms give under proper conditions contributions pro- 
portional to N*#. Neglecting these we have 


A ® «ld 2 Pp Uap + 3 2 [depp + (a mw, By)] Uap Upy + se}, (3.3) 
a, apuy 


where the presence of («u,fv) defined in (26.33) may be verified from the 
relation (26.28) between u,(k) defined in (26.23) and u,, defined in (26.21). (The 
equations refer to the book by Born and Huang.) From comparison with the 
elastic energy density (4.2) we have 


Sup = T zp, (3.4) 
Sau, pv = Tarp + (% M, BY). (3.5) 

From (2.24) we may solve d,,,g in terms of [«f, 7], 
F awe =[%B, wy] +[ar, Bu]—[a pm, By] + bap Tp» — Sev T py. (3.6) 


The result is easily checked by aid of (2.25). 
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From (2.12) and (3.4) we have 


Suy = deuve. ne dian ay. (3.7) 


Thus the complete stress tensor is expressed in force constants. We have not how- 


k | since the force constants associated 
with the boundary are involved in a basic way. 

The relation (3.7) can be useful if one makes assumptions about the force 
constants in order to obtain some model of the solid. Suppose, for example, that 
a lattice with cubic symmetry has the a priori reasonable property [6] 


LEN: LLNS 
a (by) 400 (pe) (FP) (3.8) 


I V, and y are constants and V is the volume of the solid. 


ever expressed S,,, exclusively in 0%, ( 


ne 
kk’ 
According to (3.7) we then have the following expression for the pressure P, 


where Azg ( 


* ie (3.9) 


Ve) OP ook is 
Pa, (a8 pede 
| ? aV ale. 


where OC; and C, are constants. The model represents a non-bound solid. This 
fact which easily follows in the present formalism was not noticed in the original 
treatment. 


If we describe even a modest strain with a series expansion of ® it may seem 
strongly divergent since the a ( : 
Using (2.18), however, we obtain 


eau postal) -apeo(tm lol 


4 ap 


i aa 10 a} et 
ee 
cy (iia) (to) (ea) (erg) t= B40 


it. l its 0 
where wu, b a = Uy (;) — Up (; and the atom ( ;) may be taken somewhere 


in the central part of the solid. This is still not quite convincing because of 


easily become thousands of lattice spacings. 


the fact a Ay i 
Or Ual 7p}: To see that the convergence properties are reasonable we 


consider a homogeneous deformation. The coefficient of C.2., Uap Upy Uy, iS 


l buy Ul reg 1 I 
mf 
epee (;) ee Pep, ( k’ a) ye (; .) a7) (is 4) =% 2% (;) Gapy,vh (a (3.11) 
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Equation (2.18) gives the identity 
2. Jupy,v A (;)=° (3.12) 


The corresponding relation in the limit of an infinite lattice is 


l 
2, Gepy, v2 (;) = 0. (3.13) 
If we now split (3.11) into 
l l 
82 Jaupy, vr (;) vu (+5 2 Yaupy,va (;) vu (k), (3.14) 


the first term gives a contribution essentially from the boundary region and 
both terms are reasonably expected to be proportional to N. 
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APPENDIX 
Linearized elastic constants 


Consider the work done on a solid in a deformation from positions # to + %, 


aw= |S Fiuds+ | 3 FewdV (4.1) 
a ae 


where FS and F, are the external surface and volume forces corresponding to 
the positions . From the change in the elastic energy density (see, e.g. Born— 
Huang, p. 241), 


OU 
u= > Sap Uap +4 Z, Soy, pr Uap Up + ++ Yap ~ 5 > (4.2) 
ap apy vB 
we have 
Sy, 
: bg dV vee (4.3) 


O Xp 


aw= {| wav = | > Seucmpas— [zs 
¥ 4 ap ‘ ap 


1 Now at the Quantum Chemistry Group, Uppsala University, Uppsala. 
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where n, are the direction cosines of the surface. The integrations in (4.1) and 
(4.3) are taken over the positions @. Since % is an arbitrary vector function we have 


FS = > Sap ne. (4.4) 
B 
0 Sap 
2 Rees 4.5 
F, 2 A (4.5) 


Thus S,g is the stress tensor. 
Corresponding to the positions ¢+% we have slightly different values of 


ee S.g and ng, 


(Fey = 2 Scamp. (4.6) 

The difference Sj,—S,, can be written 
Sip —Sap= > Cag. pe Upy + aisg (4.7) 

pw 
where Op, ww = Sap, ur + Sev Ou — Sap Suv (4.8) 


may be called linearized elastic constants. 
We may also use a stress tensor Szg such that 


y Sap ngd S’ = >. Sap Ng ds. (4.9) 
B B 
The difference Sip —S., can be written 


Sap — Sap = 2 S oh io Wie tee (4.10) 


Hence the coefficients Sz, in the elastic energy density may also be called 
linearized elastic constants. 


These relations may also be obtained from the theory of finite strain as 


developed by Murnaghan [7] but he does not give the linearized expressions 
explicitly. 
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